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Abstract 

The emission spectrum of a dipole embedded in a semi-infinite photonic 

crystal is calculated. For simplicity we study the case in which the dielectric 

function is sinusoidally modulated only along the direction perpendicular to 

the boundary surface plane. In addition to oscillations of the emission rate 

with the distance of the dipole from the interface we also observed that the 

shape of the emission spectrum srongly depends on the initial phase of the 

dielectric modulation. When the direction of light propagation inside the 

periodic structure is not normal to the boundary surface plane we observed 

aditional singularities in the emission spectrum, which arise due to different 

angle-dependence of the Bragg stop-band for TE and TM polarizations. 
PACS numbers: 42.70.Qs, 42.25.Gy 
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I. INTRODUCTION 



It is well known that fluorescence lifiteme of an atom can be drastically changed when 
placed in an inhomogenous medium. In his pioneering work Purcell0 predicted a strong 
enhancement of the radiative decay rate of an emitter placed inside a resonant microcavity. 
In contrast, it was also suggested! that spontaneous emission can be totally inhibited if the 
emitter transition frequency lies below that of the fundamental resonator mode. This effect 
can be understood in terms of redistribution of photonic density of states (DOS) caused 
by inhomogeneity and/or nontrivial boundary conditions imposed on the radiative field. 
During the years the spontaneous emission of a dipole coupled to various optical enviroments 
such as metallic cavitie^l, Fabry Perot two-mirror cavitiesi, dielectric microsphere^! and 
nanobubbles,i has been a subject of extensive theoretical and experimental studies.! 

Recently there has been a growing interest in studies of radiative properties of fluorescent 
molecules inside periodic dielectric structures, so called photonic crystals (PC).i'0 The Bragg 
diffraction of light that occurs in these systems opens up a spectral ga p§ (or a pseudo- 
gap) in the photonic DOS in analogy with electronic energy gaps in semiconductors. A 
quantum electro-dynamical approach for the radiative decay inside PC has revealed novel 
physical phenomena such as strong suppression of spontaneous emission within the gap 
and formation of photon-atom bound states,llHil localization of superradiant modes near 
the band edges,lii etc.. Calculations of emission spectra within the framework of classical 
theory were performed for one-dimensional Kronig-Penny type mo del0 as well as for three- 
dimensional fee lattice structuresEl. It has been established that inhibition of spontaneous 
emission within the gap is accompanied with strong enhancement at the band edges. §0 It 
was also shown that the emission spectrum strongly depends on the position of the emiter 
within the unit cell,llZllll as well as on its orientation.lll Experimental observations of inhibited 
spontaneous emission have been reported for different periodic structures. 

So far, the exsisting theoretical studies have considered infinite periodic structures. As a 
result, the emission power was identically zero within the gap.0@ In the experiments,EHH3 
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however, the photonic gap appears as a drop by a factor of ~ 2 in the emission power within 
a certain spectral interval. This points at the important role of the boundary between 
the photonic crystal and the air, which is studied in the present paper. As we will show 
below, accounting for a nearby interface results in a nontrivial dependence of the emission 
spectrum on the initial phase of the dielectric modulation. If the distance of the dipole from 
the boundary surface plane and the dielectric modulation period are respectively d and a, 
then one might expect that the dependence of the emission spectrum on the initial phase be 
small in the parameter a/d. On the contrary, we found that the strong dependence of the 
emission spectrum on the initial phase persists even in the limit d/a —>■ oo (provided there 
is no absorption in the system) . We illustrate this effect in the frame of the simplest possible 
model. Namely, we choose the dielectric modulation to be i) weak, ii) one-dimensional. Hi) 
sinusoidal. To quantitavely study the effect of a plane boundary we generalize the standard 
calculations of the emission rate in periodic media for the case of semi-infinite geometry 
(Sec. 2). In Sec. 3 we present numerical results for emission spectra illustrating the role of 
the initial phase. Discussion of our results and their relevance to recent measurements is 
presented in Sec. 4. 



II. DERIVATION OF THE POWER EMISSION SPECTRUM 

We schematically depict the system under consideration in Fig. 1. The dielectric function 
for the left half-space is constant and equals Eq, whereas for x > is given by 

e{x) — El + 5ecos{ax + (f)) , (1) 

where is the initial phase of the dielectric modulation, cr = 27r/a is the modulation wave 
vector, and 5£ is the amplitude of the modulation. Below we assume 5e <^ £i. The wave 
equations for the elctric and magnetic fields are 

V2E(r) - V(VE(r)) + ^£(a:)E(r) = t^3{r) (2) 
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V^Bfr) - V\\ne{x)] x V x B(r) + %-e{x)B{r) = — x J(r) (3) 

c 

where the radiation source J(r) = j5(r — ro) is located at point ro = {d, 0, 0). We note that 
the term proportional to 6e has been neglected in the r.h.s. of Eq. (Q). The time averaged 
radiative power per unit solid angle is given by 
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r^n{E X B*) 
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where B* is the complex conjugate of B, r = |r| and n = r/r is the unit radius vector. 
Without any loss of generality we choose = (see Fig. 1). Then it is very convenient to 
separately treat two possible orientations of the dipole. Indeed, one can easily check that the 
current density components Jz and Jx, Jy give rise to Electro-Magnetic {EM) radiation with 
respectively electric {TE polarization) and magnetic (TM polarization) fields polarized in 
the z direction. Since these two modes do not interfere, their contributions to the radiation 
power are additive. The corresponding EM wave equations for the two polarizations are 
obtained from the 2;-components of Eqs. @ and (^ by taking the Fourier transform with 
respect to y and z coordinates: 

+ - k'yjEz = -^3z5{x - d) , {TE) (5) 

where ky is the y-component of the wave vector {Ez{x; ky), Bz{x; ky) oc e*^^^). Since we want 
to calculate the power emmited in xy plane we have set A;^ = in Eqs. (|]) and (^). The 
solution of the corresponding homogenous equations may be written as a sum of incident, 
reflected and transmitted EM waves, with two linearly independent terms -^1(2;), E2{x) and 
Bi{x), B2{x) corresponding to the incident EM wave from the right and left, respectively 
(see inset of Fig. 1). To solve Eqs. (^§) we employ the variation of a constant method. 
We seek solution in the form: 

Ez{x) = Cf{x)E^{x) + C^{x)E2{x) , Bz{x) = C^{x)B^{x) + C^{x)B2{x) (7) 



Upon substituting into (|§-§D we find for tlie variational coefficients 

4:71. 
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Cf,{x) = iu—We^ / dx'E2,i{x')jJ{x' - d) (8) 



CU^) = -^^b' dx'B,^,{x') {J^ - ikyj)j5{x' - d) (9) 

where We-, Wb are the Wronskians 

dE2 dEi dB2 dBi , , 

dx dx dx dx 

and Xf'g, the lower integration limits are constants of integration. They are de- 

termined from the boundary conditions that there are no incoming EM waves, implying 
Xf = Xi = oo, X2 = X2 = —00. Hence, the solutions of Eqs. for large negative x 

which satisfy the boundary conditions can be written as follows 

47r 

E,{x) = iu—WE^3,E^{d)E2{x) , (11) 
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+iky],Bi{d))B2{x) . (12) 
i / 



c ^ V' dx 

It can be easily shown that for negative x one has Wb^E2{x) = Wb^B2{x) = ie'^^''^ /2kx- 
Thus, the remaining task is to find Ei{d), Bi{d) and {dBi/ dx)x=d- 

The solution of the homogeneous equations for the left half-space can be written as a 
sum of two plane waves, 

E^{x) = e-'^-^ + REe'^-'' , Bi{x) = e~'^-'' + RBe'^^"' (13) 

where Re-, Rb are the optical reflection coefficients for TE and TM polarizations, respec- 
tively. For the right half-space one may use the Bloch theorem to find the solution 

00 00 
E^{x) = e^«^^ ^fe^""^ , = e^"^^^ ^ e^™^ . (14) 

n=— 00 n=—oo 

When substituting Eq.(|I^) into Eqs.(^^) we obtain two infinite systems of linear, homoge- 
nous equations for the coefficents A^, 



""^1 -kl- {qE + nar)A^ + ( e^^A^, + e^^^, ) = (15) 



+ ^^"^^(^^1 - ^(^^ + + l)^)j^n+i = . (16) 

We note that the initial phase (f) exphcitly enters into these equations. 

Near the Bragg resonance that occurs at qe^b ~ o"/2, the main coefficents which con- 
tribute to the sums in Eq.(|lJ) are , Aq and A^^, A^^ since the rest of A^'^ are small in 
the parameter 6e/ei. In this approximation the equation systems (p!^- [l6|) may be simplified 
into two 2x2 matrix equations. Requiring for the determinants of these matrices to vanish, 
one finds the dispersion relations for the two EM polarizations near resonance 
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where we have introduced A = fe/4ei, qe^b — o"/2 = SqE,B ^ cr/2, u — ujq = Su <ti ujq, 
K = a^c^/AujQei and the fcy-dependent resonant frequency c^o is given by 



a;o = -V^V4 + A;2. (19) 



Eqs. (p!7|-|T8D show that there is a spectral gap for EM waves propagating in the system 
centered at Uq. For both TE and TM polarizations the central gap position shifts to 
higher frequencies with increasing ky (which also determines the propagation direction of 
the radiative field). For TE polarization the gap broadens with increasing angle 6' (see 
Fig. 1) whereas for TM polarization the gap narrows and disappears at 2k = 1, which 
corressponds to the propagation direction for which ky = a/2. This can be defined as a 
Brewster angle for Bragg diffraction. If one increases ky further then the gap reopens again. 

Using Eqs. (P^-(P!^) we obtain the following expressions for the electric and magnetic 
fields for the right half space: 



Ei{x) = A^^e'^^^^ ( e^?^ - FEe-''''e-'^^] (20) 



5i(x) = A^ie*^^^^ ( e*5^ - F^e-^^e-^^^ ) (21) 



Here the functions Fe and Fb describe the couphng between incident and Bragg reflected 
waves and are defined as follows: 



Fs(S^,e) = —^C-^-JC-^)'-AHl-2nr]. (23) 



A(l — 2k) Kluq \I Kluq, 

Finally, we match the fields components and their derivatives at x = 0, take the inverse 
Fourier transform with respect to ky,kz and use Eq. (^ to get the following expression 
for the total radiated power dP/dQ, normalized to the radiation power when there is no 
dielectric modulation: 

'dP\ 



- 



Jz-'-OE 


1 _ p^Q-i(^'i+4,) 




1 - RoeFec-^'^ 



1 - xFBe-M+<^) ^,s,,d 



-e 



(24) 



1 - RoEFse-'^ 

Here Tqe^ Roe and Tqb, Rqb are the conventional Fresnel's transmission and reflection 
coefficients from the dielectric interface without the dielectric modulation for the amplitudes 
of TE and TM polarized waves, respectively, ji = jy cos 9' + jx sin 6', where 6' is related 
to the observation angle 6 by Snell's law. The quantity x? which depends on the dipole 
orientation in the xy plane is defined as x = {jy — jx tan 9') / {jy + j^ tan 9'). 

Eq. (p^ ) is the main result of this paper. We note that the phase is explicitly present in 
both the numerator and denominator of Eq. (p^, indicating the important role of the initial 
modulation phase at the boundary interface. In the next section we numerically analyze the 
role of (f) in the emission spectrum for different situations. 
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III. NUMERICAL RESULTS 



Before presenting our numerical results let us concentrate on a particular experimental 
case of opalsH'il and opal replica.il These artificial photonic crystals consist of closely packed 
Si02 spheres forming an fee structure that contains fully interconnected voids. The opal 
replica is formed by filling these voids with a precursor polymer solution and then etching 
the 5*^02 spheres after polymerization.ii The opals and opal replica PC are infiltrated with 
various fluorescent dye solutions to provide a radiation source inside the crystal. Inhibited 
spontaneous emission of the dye molecules in such PC have been recently studied by several 
groups.ii'ii The refractive index n of Si02 is n ~ 1.46 and the refractive index contrast An 
between the Si02 balls and the dye solution ranges between An ^ 0.1 — 0.3. This is not 
sufficient for a formation of a complete photonic band gap. Instead, the system posseses 
pseudogaps (or partial gaps) with an angle-dependent central frequency. To compare our 
results to the experiments we have chosen in our model A = 5e/Aei = 0.1. 

First, we consider the case when the emitter is many periods away from the interface 
[d = 5a). 1^ In Fig. 2 we show the emission power averaged over the orientations of the 
emitter as well as over its position within the unit cell (the zero on the frequency axis in 
all plots corresponds to the central gap frequency Qq at 9 = 0, where Qq = (7c/{2^yei)). 
Inside the gap, the emission power is strongly suppressed. It is seen, however, that even in 
the limit of large d/a the features of the averaged emission power still depend on the initial 
phase. 

The evolution of {dP/dQ)^ with (p outside the gap is described as follows. At = 
there is a well pronounced singularity at the lower edge of the spectral gap. With increasing 
(f), this singularity diminishes whereas another singularity starts to develop at the upper 
edge; the spectrum becomes symmetric at = 7r/2. Further increase of leads to a gradual 
transformation of the initial curve into its mirror image with respect to the central gap 
frequency, the singularity now occuring at the upper band edge. 

Fig. 2 corresponds to 6* = 0. The sensitivity of the emission power to (p appears to be 
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even more pronounced for ^ 7^ 0. We illustrate this effect by plotting in Fig. 3 the averaged 
{dP/(Kl)j^ for 9 = 60" for the set of initial phases = 0, 7r/4, 7r/2. With increasing (p again, 
it is seen that there is a tendency for the emission spectrum to become symmetric near 
(j) = 71/2. At0 = 7r (not presented here) the emission spectrum is again transformed into 
the mirror image of the initial curve at = but now with respect to the shifted central 
gap frequency VLq/cos9'. We also note the appearance of additional singularities in the 
emission spectrum. Their origin lies in the different angle-dependenceies of the gap width 
for TE and TM polarizations, as seen before in Eqs (0) and ([Tsl). To be more specific, 
we note that the highest and lowest frequency peaks correspond to the band edges for TE 
polarization. Similarly, two peaks at the intermidiate frequencies determine the band edges 
for TM polarization. Both pairs of the singularities are located symetrically around the 
shifted Bragg frequency ujq = VLq/ cos 9'. 

Let us now turn to the discussion of the case where the emitter is close to the interface. 
The main feature of this situation is that {dP/ dQ) ^ inside the gap remains finite. One can 
see from Eq. ( [2^ ) that moving the emitter N periods away from the boundary decreases the 
emission power at the center of the gap by a factor of exp{—2iTNA). Therefore, to study the 
features of the spontaneous emission inside the gap we choose N = 1,2. This situation is 
illustrated in Fig. 4 where we plot the averaged {dP/ dQ) n for the observation angle 9 = 60°. 
One can see that the evolution of the emission power with increasing the initial phase is 
very strong. Again, we note that there are four well pronounced singularities in the emission 
spectrum that corresponds to the spectral gap edges of the TE and TM polarizations. In 
particular, for = there is a noticable enhancement of the emission rate (by a factor 
of ~ 2) at the frequency which determines the lower band-edge of the TE polarization, 
whereas for = 7r/2 a similar enhancement occurs at the edges of the TM polarization gap. 

Finally, in Fig. 5 we plot the emission power integrated over the observation angle 9 (i.e., 
the total power emitter in xy plane) for = and = vr. In this case too an averaging over 
the dipole position within the unit cell has been performed. To allow for an unpolarized 
emission we have chosen jx = jy = jz- Remarkably, even after angular averaging a weak 
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dependence of the emission power on the initial phase still persists. 

IV. DISCUSSION 

We have studied the emission spectrum of a dipole inside a one-dimensional periodic 
structure in the presence of a nearby plane boundary. As expected, the emission rates are 
strongly suppressed for frequencies inside the spectral gap, provided that the emitter is many 
periods away from the interface {d/a > 5). For frequencies near the band edges, the emission 
spectrum changes drastically with the initial phase of the dielectric modulation. We also 
observed enhancement of the emission rates near the band edges, however, by a factor much 
smaller than predicted in the previous studies.0lli This can be attributed to the fact that 
the modulation in our model is weak. In Ref. where the total radiation power from a 
dipole inside a infinite 3D fee lattice was calculated numerically, the contrast in the dielectric 
constant Ae was more than 10 allowing for the formation of a complete photonic band gap. 
This resulted in an enhancement of the radiated power at the band edges by a factor of ~ 25. 
Similarly, in Ref. |1^ where the authors considered ID Kronig-Penney type modulation of the 
refractive index, the enhancement factor at the band edges was about ~ 30, whereas inside 
the gap it was identically zero. Apparently, this resulted from consideration of radiative 
modes polarized parallel to the dipole direction. Allowing for nonpolarized radiation [i.e. in 
all directions) should lead to qualitatively different results (see for example Fig. 5). Strictly 
speaking, the straightforward comparision of our calculation with those of previous studies 
is not possible due to the different approach developed here. However, it is clear from our 
consideration that the sensitivity of the power emission to the boundary conditions should 
persist also for strong and non-sinusoidal modulation and become even stronger. 

In this paper we also studied the case where the emitter is sufficiently close to the 
interface, so that the emission power for the frequencies inside the gap is finite. In this case 
also we showed that the features of the emission spectrum are very sensitive to the initial 
phase of the periodic modulation, hence emphasizing the effect of the boundary. 
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We note that in our calculations we assumed that no deffects exist in the system. In 
the presence of weak disorder the emission spectrum would be significantly modified. Let 
us consider radiation from a dipole that is many periods away from the boundary with 
frequency inside the spectral gap. If there are no defects then the radiation (e. g., in the 
direction normal to the boundary) is strongly attenuated (see Fig. 2). However, introducing 
a small concentration of defects opens up a new mechanism for the light emission to come 
out from the PC in the direction normal to the boundary. Namely, light can propagate 
without any attenuation in directions for which the Bragg condition is not satisfied and 
then scatter off defects that are close to the interface. As a result, the emission power in the 
direction normal to the boundary is not exponentially small but is finite, similar to the case 
of a dipole close to the boundary. Remarkably, only defects close to the boundary, within 
the Bragg atenuation lenght = (27rcrA)~^ from the interface, contribute to the emission 
power. Owing to this effect, we conclude that the emission spectrum in the presence of a 
weak disorder is rather universal, since it does not depend on the dipole-interface distance d. 
Experimentally this was demonstrated in Ref. ^ where the authors found a way to excite 
fluorescent dye molecules at different distances from the boundary. They have demonstrated 
that the emission power at the center of the spectral gap is attenuated by a factor of ~ 2 
and does not change upon changing d. This was interpretated in terms of light scattering 
off defects close to the PC interface. 
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FIGURES 

FIG. 1. Schematic representation of our PC system. The emitter is located on the x axis at a 
distance d from the interface. The miit vector n hes in xy plane. Inset: Illustration of two linearly 
independent solutions of the homogenous equations discussed in the text. 

FIG. 2. Normalized emission power {dP/dQ)^ averaged over the dipole emitter orientations at 
= for ^ = 0, 7r/2 and tt. Here 5uj = u — VLqis the deviation from the Bragg frequency at = 0. 
Calculations were performed with ei = Seq = 3 and A = 0.1. The emitter is five periods away 
from the interface {d = 5a) implying strong suppression of the emission within the gap. 

FIG. 3. Normalized and averaged power emission {dP/dVL)j^ at 9 = 60" for (p = 0, 7r/4 and 7r/2. 
The rest of the parameters are the same as in Fig. 2. 

FIG. 4. Emission spectrum of the emitter located at (a) iV = 1, (b) AT = 2 periods away from 
the boundary for 9 = 60° and for the same set of the initial phases as in Fig. 3. Averaging over 
the position within the unit cell and over all possible orientaions of the dipole emitter direction 
has been performed. 

FIG. 5. The total power emitted in xy plane for cf) = (solid line), (f) = it (dashed line). Here 
3x = jy = jzi d = 2a. Averaging over the dipole emitter position within the unit cell has been 
performed. 
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